Acoustic phonon scattering limited carrier mobility in 2D extrinsic graphene 
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We theoretically calculate the phonon scattering limited electron mobility in extrinsic (i.e. gated 
or doped with a tunable and finite carrier density) 2D graphene layers as a function of temperature 
(T) and carrier density (n). We find a temperature dependent phonon- limited resistivity p P h{T) 
to be linear in temperature for T > 50K with the room temperature intrinsic mobility reaching 
values above 10 5 cm 2 /Vs. We comment on the low-temperature Bloch-Griineisen behavior where 
Pph(T) ~ T 4 for unscreened electron-phonon coupling. 

PACS numbers: 81.05.Uw, 72.10.-d, 73.40.-c 



I. INTRODUCTION 

Low temperature carrier transport properties of 2D 
graphene layers have been of great current interest to 
both experimentalists^^^ and theorists^ ' 7 i 8 1 9 1 10 i alike 
ever since the possibility of fabricating stable gated 2D 
graphene monolayers on SiC>2 substrates and measur- 
ing the density-dependent conductivity of the 2D chi- 
ral graphene carriers were demonstrated!'. Much of the 
early interest focused understandably on the important 
issues of the scattering mechanisms limiting the low- 
temperature conductivity and the associated graphene 
"minimal conductivity" at the charge neutral ("Dirac") 
point. One of the dominant low-temperature scattering 
mechanisms^^iiS in graphene is that due to screened 
Coulomb scattering by unintended charged impurities in- 
variably present in the (mostly SiC>2) substrate (and the 
substrate-graphene interface) although short-range scat- 
tering by neutral defects also contributes, particularly at 
high carrier densities in high-mobility samples. It has 
therefore been argued*^ that gated graphene layers are 
similar to 2D electron systems in confined semiconductor 
structures (e.g. Si inversion layers, GaAs heterostruc- 
tures and quantum wells) where also long-range charged 
impurity scattering dominates low-temperature ohmic 
transport with short-range (e.g. interface roughness) 
scattering playing a role at high carrier densities! 11 ' 12 

Given the exciting technological context of graphene as 
a prospective electronic transistor material for future ap- 
plications, the question therefore naturally arises about 
the limiting value of the intrinsic room-temperature 
graphene mobility if all extrinsic scattering mechanisms, 
e.g. charged impurities, neutral defects, interface rough- 
ness, graphene ripples etc., can be eliminated from the 
system. This question is more than of academic interest 
since serious experimental efforts are underway^ to elim- 
inate charged impurities from graphene by using differ- 
ent substrates or by working with free standing graphene 
layers without any substrates. It is also noteworthy that 
the systematic elimination of charged impurity scatter- 
ing through modulation doping and materials improve- 
ment in the MBE growth technique has led to an as- 
tonishing 3,000-fold enhancement in the low-temperature 



2D GaAs electron mobility from 10 4 cm 2 /Vs in 1978 to 
30 x 10 6 cm 2 /Vs in 2000, future enhancement to 100 
million cm 2 /Vs mobility is anticipated^ in the next few 
years. 



One great advantage of graphene over high-mobility 
2D GaAs systems is that the lack of strong long- 
range polar optical phonon scattering, which completely 
dominates 1 ^ the room-temperature GaAs mobility (~ 
2, 000 cm 2 /Vs), in graphene should lead to very high in- 
trinsic room temperature graphene mobility, limited only 
by the weak deformation-potential scattering from the 
thermal lattice acoustic phonons. In this work, we cal- 
culate the temperature-dependent 2D graphene mobility 
limited only by the background lattice acoustic phonon 
scattering. We find that room-temperature intrinsic (i.e. 
just phonon-limited) graphene mobility surpassing 10 5 
cm 2 /Vs is feasible using the generally accepted values 
in the literature for the graphene sound velocity and de- 
formation coupling. There is some uncertainty in the 
precise value of the electron-phonon deformation poten- 
tial coupling constant, leading to a concomitant uncer- 
tainty in the intrinsic graphene mobility. This situation 
is similar— to the 2D GaAs system where in fact pre- 
cise measurement of the phonon-limited 2D mobility led 
to the correct deformation potential coupling for trans- 
port studies, and this could be the case for graphene 
also where a quantitative comparison between our theo- 
retical results presented in this work with the measured 
temperature-dependent graphene mobility, very recently 
becoming availabl o 17 ' 18 , could lead to an accurate deter- 
mination of the graphene electron-phonon deformation 
potential coupling constant. 



The paper is organized as follows. In section II the 
Boltzmann transport theory is presented to calculate 
acoustic phonon scattering limited 2D graphene conduc- 
tivity. Section III presents the results of the calculation. 
In section IV we discuss the results compared to experi- 
mental data, and we conclude in Section V. 
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II. THEORY 

We use the Boltzmann transport theor y 15 i 16 to cal- 
culate acoustic phonon scattering limited 2D graphene 
conductivity. We consider only the longitudinal acoustic 
(LA) phonons in our theory since either the couplings to 
other graphene lattice phonon modes are too weak or the 
energy scales of these (optical) phonon modes are far too 
high for them to provide an effective scattering channel 
in the temperature range (5 — 50QK) of our interest. 

The conductivity of graphene is given by 



e 2 D(E F )f(r) 



(1) 



where vp is the Fermi velocity, D(Ep) = 
(gs9v/^fi 2 )EF/vp is the density of states of graphene 
at the Fermi level (Ep) and (r) is the relaxation time 
averaged over energy, i.e. 



(2) 



where /(e) is the Fermi distribution function, f(ek) — 
{1 + cxp[/3(e fc - n)}}- 1 with (3 = l/k B T and ju(T,n) 
as the finite temperature chemical potential determined 
self-consistently. The energy dependent relaxation time 
[r(ek)] is defined by 
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where #kk' is the scattering angle between k and k', 
e = Hvp |k|, and Wkk' is the transition probability from 
the state with momentum k to k' state. In this paper 
we only consider the relaxation time due to deformation 
potential (DP) coupled acoustic phonon mode. The de- 
formation potential due to quasi-static deformation of 
lattice is taken into account. Then the transition proba- 
bility has the form 



2tt 



W^ = T ^2\C(q)\ 2 A(e > e') 



(4) 



where C(q) is the matrix element for scattering by acous- 
tic phonon and A(e,e') is given by 

A(e, e') = N q 5(e -e' + w q ) + (N q + l)S(e -s'- w q ), (5) 

where a-> q = v p h(\ is the acoustic phonon energy with v p h 
being the phonon velocity and N q is the phonon occupa- 
tion number 



9 exp(/3w q ) - 1 ' 



(6) 



The first (second) term is Eq. corresponds to the ab- 
sorption (emission) of an acoustic phonon of wave vector 



q = k — k'. The matrix element C(q) is independent 
of the phonon occupation numbers. The matrix element 
|C(q)| 2 for the deformation potential is given by 



l^(q)l 2 



D 2 hq 
2Ap m v p h 



1 - 



\2k) 



(7) 



where D is the deformation potential coupling constant, 
p m is the graphene mass density, and A is the area of the 
sample. 

The scattering of electrons by acoustic phonons may be 
considered quasi-elastic since hujq <C Ep, where Ep is the 
Fermi energy. There are two transport regimes, which ap- 
ply to the temperature regimes T <C Tbg and T 3> T^g, 
depending on whether the phonon system is degener- 
ate (Bloch-Griineisen, BG) or non-degenerate (equipar- 
tition, EP). The characteristic temperature Tbg is de- 
fined as ksTBG = 2kpv p h, which is given, in graphene, 
by Tbg — ^Vphkp/kB ~ 54y / n K with density mea- 
sured in unit of n = 10 12 cmT 2 . First we consider 
Su q <C ksT. In this case we have N q ~ ksT '/huj q , and 
A(e,e') = (2k B T/Hu}q)S{e - e'). Then the relaxation 
time is calculated to be 



1 e k D 2 



r(e k ) h 3 Av 2 F p„ 



-k B T. 



(8) 



J P h 



Thus, in the non-degenerate EP regime (fiw q -C ksT) the 
scattering rate [l/r(ek)] depends linearly on the temper- 
ature. Since at low temperatures (Tbg *C T Epjks) 
(t) ps t (Ef) the calculated conductivity is indepen- 
dent of Fermi energy or equivalently the electron density. 
Therefore the electronic mobility in graphene is inversely 
proportional to the carrier density, i.e. p cx 1/n. The EP 
regime has recently been considered in the literatur e 18 ' 19 . 
We note that the similar linear temperature dependence 
of the scattering time has been reported for nanotubes^ 
and graphites 21 . 

To calculate the relaxation times in the BG regime 
where hujq ~ ksT we have to keep the full form as in Eq. 
©. Since the acoustic-phonon energy is comparable to 
fcflT the temperature dependence of the relaxation time 
via the statistical occupation factors in Eq. © becomes 
more complicated. In BG regime the scattering rate is 
strongly reduced by the occupation factors because for 
phonon absorption the phonon population decreases ex- 
ponentially and also phonon emission is prohibited by 
a sharp Fermi distribution. To calculate the low tem- 
perature behavior of the resistivity we can rewrite the 
averaged inverse scattering time over energy as 



1 2E F 



dq(l-cos6)\C(q)\ 2 G(u q ), (9) 



2TTh(hv F ) 2 1 ™"' rwi ~ v q 

where q = 2kp s'm(0/2) and G(oj) is given by22, 
1 



k B T 
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def(s){N q [l-f(e + u>)] 

+(JV 9 + l)[l-/(£-w)]} 
N q (N q + 1). (10) 
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FIG. 1: (Color online) Calculated inverse relaxation times as 
a function of energy for different temperatures T/Tbg = 0.2, 
0.5, 1.0, and 1.5 for an electron density n = 10 12 cm" 2 with 
Tbg = 54ff. The deformation potential coupling constant 
D = 19 eV and the phonon velocity v p h — 2 x 10 6 cm/s are 
used in this calculation. 



Then we have in low temperature limits T <C Tbg 



1 



1 1 



D 2 



4!C(4) 
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Thus, we find that the temperature dependent resis- 
tivity in BG regime becomes p ~ T 4 without screen- 
ing effects. If we include screening effects by the carri- 
ers themselves 2 ^ the low-temperature resistivity goes as 
p ~ T 6 . The screening effects on the bare scattering rates 
can be introduced by dividing the matrix elements C (q) 
by the dielectric function of graphene. But the matrix el- 
ements in graphene arise from the change in the overlap 
between orbitals placed on different atoms and not from 
a Coulomb potential. Thus, we neglect screening effects 
in the calculation, and only consider unscreened defor- 
mation potential coupling. Even though the resistivity 
in EP regime is density independent, Eq. indicates 
that the calculated resistivity in BG regime is inversely 
proportional to the density, i.e. pbg ~ 1/™, or cquiva- 
lently the mobility in BG regime is density independent. 



III. RESULTS 

In this calculation we use the following parameters: 
graphene mass density p m = 7.6 x 1CT 8 g/cm 2 , acoustic 
phonon velocity v p h = 2 x 10 6 cto/s, and deformation po- 
tential D = 19eV. Even though the phonon velocity v p h 
is well defined experimentally the value of the deforma- 
tion potential coupling constant is not establishe d 21 ' 24 . 
In general, the constant D could be obtained on the ba- 
sis of the fact that the shift of energy dispersion from its 
equilibrium state reaches the order of the atomic energy, 
i.e. D ~ e 2 /a with a being the lattice constant, which 
is of the order of 10 eV in graphene. We note that we 
have used ref. 24 for obtaining the phonon parameters, 
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FIG. 2: (Color online) (a) Calculated graphene resistivity as a 
function of temperature for several densities n = 1, 3, 5xl0 12 
cm -2 . We use the deformation potential D = 19 eV. Note 
that in BG regime (T < Tbg) P ~ T 4 and in EP regime 
(T > WOK) p ~ T. (b) The same as Fig. 2(a) in linear 
scale. Inset shows the resistivity in low temperature limits 
(T < 80K). 



but different values of D, differing by factors of three 
(i.e. D w 10 — 30 eV), are quote d 19 ' 21 ' 24 in the literature. 
Since r _1 oc D 2 , the resulting graphene resistivity could 
differ by an order of magnitude depending on the precise 
value of D. 

In Fig. 1 we show the calculated inverse relaxation 
times for deformation potential scattering by acoustic 
phonon as a function of energy for different temperatures 
T/Tbg = 0.2, 0.5, 1.0, and 1.5 for an electron density 
n = 10 12 cm" 2 with Tbg — 54JT. The inverse relaxation 
time in BG regime (T < Tbg) shows a characteristic dip 
(suppression of scattering rate) in a narrow region around 
Fermi energy Ep due to the statistical occupation factors. 
Above Bloch-Griineisen temperature (T > Tbg) the dip 
structure disappear and the scattering rate becomes close 
to the scattering rate of equipartition regime. 

In Fig. 2(a) we show our calculated graphene resistiv- 
ity, p = a~ l , as a function of temperature on a log-log 
plots, clearly demonstrating the two different regimes: 
BG p - T 4 behavior for T < T BG ~ 20 - WOK and the 
EP p - T behavior for T > lOOif. We note that T BG 
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FIG. 3: (Color online) Calculated mobility limited by the 
acoustic phonon with the deformation potential coupling con- 
stant D = 19 eV (a) as a function of temperature for different 
densities n = 1, 3, 5xl0 12 cm -2 and (b) as a function of den- 
sity for different temperatures T=77K and T=300K. 



(oc y/n) depends weakly on density, and the true BG be- 
havior is likely to show up at relatively low temperatures. 
In Fig. 2(b), we show p{T) for several densities on a lin- 
ear plot, emphasizing the strong linear in T dependence 
of the acoustic phonon-limited graphene resistivity rang- 
ing from 100K to 500K. This rather large temperature 
range of p(T) ~ T behavior of acoustic phonon scattering 
limited resistivity is quite generic to 2D semiconductor 
structures, and our finding for graphene here is qualita- 
tively similar to what was earlier found to be the ease 
for 2D GaAs structure a 15 ' 16 . The crucial difference be- 
tween graphene and 2D GaAs is that in the latter system 
polar optical phonon scattering becomes exponentially 
more important for T > 100A and dominates at room 
temperatures whereas in graphene we predict a linear 
2D resistivity upto very high temperatures (~ 1000 A") 
since the relevant optical phonon have very high energy 
(~ 2000A) and are simply irrelevant for carrier trans- 
port. 

In Fig. 3, we show our calculated intrinsic graphene 
mobility, p = (enp) , as functions of temperature and 
carrier density. Within our model, the unscreened acous- 
tic phonon scattering limited graphene mobility is in- 



versely proportion to T and n individually for T > 100A. 
Assuming D = 19eV, as used in Fig. 3, p could reach 
values as high as 10 5 cm 2 /Vs for lower carrier densities 
(n < 10 12 cm -2 ). For larger (smaller) values of D, p 
would be smaller (larger) by a factor of D 2 . It may be 
important to emphasize here that we know of no other 
system where the intrinsic room-temperature carrier mo- 
bility could reach a value as high as 10 5 cm 2 /Vs. 



IV. DISCUSSION 

The three key theoretical findings on phonon-limited 
graphene mobility of this work are : (1) p ~ T for 
T > 100 A; (2) p cx T 4 (T 6 ) for T < 50 A) for un- 
screened (screened) deformation potential coupling; (3) 
I 1 3.7 x 10 7 /D 2 h cm 2 /Vs at room temperature (T = 
300A) where D is measured in eV and n is carrier den- 
sity n measured in units of 10 12 cm~ 2 . These theoretical 
predictions being rather precise, the question naturally 
arises about the experimental status and the verifica- 
tion of our theory. Very recently, experimental graphene 
transport data at room temperature (or even abovei 7 -) 
have started becoming availableA 7 -^. The only aspect 
of our theory that can be directly compared with the 
existing experiment is the p cx T behavior at high tem- 
peratures (> 200 A), and this is indeed consistent with 
the recent data from two different group a 17 i 18 . Geim and 
collaborators have recently concluded-^ that the room 
temperature graphene intrinsic mobility could be as high 
as 10 5 cm 2 /Vs, which is also consistent with our the- 
ory. But, as emphasized by us, the actual mobility value 
varies inversely as D 2 , and therefore a precise knowledge 
of the deformation potential coupling is required for an 
accurate estimate of the intrinsic mobility. 

A more detailed comparison between our theory and 
the experimental results on p(T) shows some qualita- 
tive difference which are not understood at this point. 
For example, the experimental crossove r 17 i 18 to the high- 
temperature linear (p ~ T) behavior in the intrinsic re- 
sistivity appears to be closer to a T 2 behavior rather 
than the T 4 BG behavior we predict. More disturbingly, 
the experimental crossover from the high-temperature 
linear behavior to the low-temperature high power law 
behavior appears to be occurring at a much higher tem- 
perature (100 — 200A) than the theoretical prediction 
(20 — 50 A). At this stage we have no explanation for 
the lower-temperature disagreement between experiment 
and theory, but below we discuss several possibilities. 

The experimentally measured resistivity in the cur- 
rent graphene samples is completely dominated by ex- 
trinsic scattering (and not by phonon scattering) even 
at room temperatures since the low-temperature (< 4A) 
mobility is typically 5, 000 — 15, 000 cm 2 /Vs, and the 
intrinsic room-temperature phonon contribution, as ob- 
tained theoretically by us or inferre d 17 ' 18 from recent 
temperature-dependent experiments, is 10 — 20 times 
larger (~ 100, 000cm 2 /Vs). This means that any ex- 
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FIG. 4: (Color online) The deviations from Mattiessen's rule 
for two different densities, (a) n = 10 12 cm -2 and (b) n — 
7 x 10 11 cm -2 . Here pi (p p h) represents the resistivity due to 
impurity (phonon) scattering. 



perimental extraction of the pure phonon contribution 
to graphene resistivity involves subtraction of two large 
resistances (i.e. the measured total resistance and the 
extrapolated T — extrinsic temperature-independent 
resistance arising from impurity and defect scattering) 
of the order of kfl each to get a phonon contribution 
roughly of the order of 100 fi. Apart from the inher- 
ent danger of large unknown errors involved in the sub- 
traction of two large numbers to obtain a much smaller 
number associated with phonon scattering contribution 
to graphene mobility, there is the additional assumption 
of the Matthiessen's rule, i.e. p tot = p p h+Pi where p tot is 
the total resistivity contributed by impurities and defects 
(pi) and phonons (p p h), which is simply not valid. In par- 
ticular, the impurity contribution to resistivity also has 
a temperature dependence arising from Fermi statistics 
and screening which, although weak, cannot be neglected 
in extracting the phonon contribution (particularly since 
the total phonon contribution itself is much smaller than 
the total extrinsic contribution). In particular, the tem- 
perature dependent part of the charged impurity scatter- 
ing contribution to graphene resistivity could be positive 



or negative^ depending on whether screening or degener- 
acy effects dominate, and therefore the phonon contribu- 
tion, as determined by a simple subtraction, could have 
large errors, particularly in the low (T < 1000/f ) tem- 
perature regime. Indeed, a recent measurement^ of p(T) 
in the — regime finds small temperature depen- 

dent contributions to graphene resistivity which could be 
either positive or negative depending on the sample mo- 
bility and which, in all likelihood, arises from extrinsic 
impurity scattering. 

In Fig. 4 we show the failure of the Matthiessen's 
rule in graphene (particularly at higher/lower temper- 
atures/densities) by calculating the total graphene re- 
sistivity arising from screened charged impurity scatter- 
ing (pi) and phonon scattering (p p h) — it is clear that 
Ptot 7^ P P h + Pi for lower/higher densities/temperatures. 

If the experimentally extracte d 17 ! 18 phonon contribu- 
tion to the graphene resistivity turns out to be accurate 
in spite of the rather questionable subtraction procedure 
discussed above, then the disagreement between our in- 
termediate temperature (50 — 100K) theoretical results 
and the experimental data would indicate the presence 
of some additional phonon modes which must be partic- 
ipating in the scattering process. We can, in fact, get 
reasonable agreement between our theory and the exper- 
imental data by arbitrarily shifting the BG temperature 
Tbg to a higher temperature around 200K. This shift 
could indicate a typical phonon scale which causes ad- 
ditional scattering other than the LA phonons coupled 
to the carriers through the deformation potential cou- 
pling considered in our work. At this stage we cannot 
speculate on what these additional modes could be. One 
possibility is that these are the zone-edge out of plane ZA 
phonon modes with vibrations transverse to the graphene 
plane 1 ' . (In the Appendix we provide the calculated car- 
rier resistivity in the presence of an additional phonon 
mode with a soft gap.) Another possibility considered in 
ref. [TH is that these are the thermal fluctuations ("rip- 
plons") of the mechanical ripples invariably present in 
graphene samples^. Of course such additional "phonon" 
scattering channels will lead to additional unknown cou- 
pling parameters making the resultant theory essentially 
a data fitting procedure. The advantage of our minimal 
theory is that it involves only two phonon parameters: D 
and v p h associated with the 2D graphene LA phonons. 
More data in higher mobility samples will be needed to 
settle this question since the subtraction problem inher- 
ent in the current technique for extracting the phonon 
contribution would make analyzing this issue a difficult 
task. 



V. CONCLUSION 

We have calculated the intrinsic temperature depen- 
dent 2D graphene transport behavior upto 500K by 
considering temperature and density dependent scatter- 
ing of carriers by acoustic phonons. We have provided 
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FIG. 5: (Color online) Calculated resistivity with both the 
acoustic phonon scattering and the inter-valley phonon scat- 
tering. See appexdix for details. 

a critical discussion of our results in light of the re- 
cent experiment s 1 7 i 18 . The lack of precise quantitative 
knowledge about graphene deformation potential cou- 
pling makes a quantitative comparison with the exper- 
imental data problematic. 

We thank Michael Fuhrer (ref. [T3|) and Andre Geim 
(ref. 18]) for sharing with us their unpublished data, 
and Andre Geim for a careful reading of our manuscript. 
This work was supported by U.S. ONR, LPS-NSA, and 
SWAN-NSF-NRI. 



APPENDIX 

In this appendix we calculate the phonon scattering 
limited carrier mobility including effects of two phonon 
branches: the regular LA phonon (as considered in the 
main part of this paper) and an additional "intervalley" 
phonon branch with a soft gap (~ 70 meV) representing 
the inter-valley phonon, the ZA phonon mode at the K 
point. The theoretical motivation is to demonstrate that 
the combination of the LA phonon and an optical phonon 
(i.e. the inter-valley ZA mode) with a softy gap could 
indeed lead to qualitative (or even quantitative, if the 
phonon parameters have appropriate values) agreement 
between theory and experiment. 

In this context we consider the long wavelength phonon 



scattering. As temperature increases phonons of large 
wave vectors are involved in the scattering in multi- 
valley structures. Thus the inter-valley phonon scat- 
tering becomes significant at high temperatures^. In 
graphene there are two minima of the conduction band 
at K and K' points in Brillouin zone. The scattering 
between K and K' points requires the participation of 
inter- valley phonons, whose wave vectors are close to 
qy = kx — kx' and the frequencies of these phonons are 
close to u>ij — Lu(qij). The relaxation time for inter- valley 
phonon scattering may be considered by assuming con- 
stant inter- valley phonon energies fay . Then the matrix 
element for the inter-valley phonon scattering becomes 
|C(q)| 2 = hD^j/2Ap m 0Jij where Djj is the deformation 
potential coupling constant for inter-valley phonons in 
unit of eV/A. Since toy <C Ep in graphene the scatter- 
ing of electrons from inter-valley phonons is considered 
quasielastically. 

In Fig. 5 we show the calculated resistivity with 
both the acoustic phonon scattering and the inter-valley 
phonon scattering as a function of temperature for two 
different densities. The following parameters are used 
in this calculation: deformation potential coupling con- 
stant D — lOeV, acoustic phonon velocity v p h = 2 x 10 6 
cm/s, Dij — 7 eV/A, and inter-valley phonon energy 
tiWij — 70 meV which corresponds to the lowest phonon 
energy (ZA) at the K point. The calculated resistivities 
have very weak density dependence, which comes from 
the energy averaging. Below 200K the acoustic phonon 
scattering dominates (linear in temperature), but above 
200 K both phonon scatterings contribute in the trans- 
port. Note that the temperature dependence of the resis- 
tivity is linear in both regimes, but has different slopes. 
The high power law behavior (p ~ T 4 ) only applies at 
very low temperatures (T < 50 K). Basically, there is 
a sharp turn-on in phonon scattering in the 150-250K 
range as the inter-valley ZA phonon scattering becomes 
effective. 

Whether the results shown in Fig. 5, including the 
effects of inter-valley phonon scattering are physically 
meaningful or not will depend on the direct observation 
of these ZA phonon modes via Raman scattering exper- 
iments. At this stage, all we have established is that in- 
clusion of this additional soft-gap phonon mode gives im- 
pressive agreement between theory and exp er imen t 1 7 1 1 8 . 
More work is needed to validate the model of combined 
LA and ZA phonon scattering contributing to the tem- 
perature dependent graphene resistivity. 
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